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ABSTRACT 
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In this work we report on high-resolution IR absorption studies that provide a detailed view 
on how the peripheral structure of irregular polycyclic aromatic hydrocarbons (PAHs) affects the 
shape and position of their 3-wm absorption band. To this purpose we present mass-selected, 
high-resolution absorption spectra of cold and isolated phenanthrene, pyrene, benz|alantracene, 
chrysene, triphenylene, and perylene molecules in the 2950-3150 cm~! range. The experimental 
spectra are compared with standard harmonic calculations, and anharmonic calculations using 
a modified version of the SPECTRO program that incorporates a Fermi resonance treatment 
utilizing intensity redistribution. We show that the 3-um region is dominated by the effects of 
anharmonicity, resulting in many more bands than would have been expected in a purely harmonic 
approximation. Importantly, we find that anharmonic spectra as calculated by SPECTRO are in 
good agreement with the experimental spectra. Together with previously reported high-resolution 
spectra of linear acenes, the present spectra provide us with an extensive dataset of spectra of 
PAHs with a varying number of aromatic rings, with geometries that range from open to highly- 
condensed structures, and featuring CH groups in all possible edge configurations. We discuss 
the astrophysical implications of the comparison of these spectra on the interpretation of the 
appearance of the aromatic infrared 3-~m band, and on features such as the two-component 
emission character of this band and the 3-jzm emission plateau. 


Subject headings: astrochemistry — ISM: molecules — methods: laboratory — techniques: spectroscopic 


— line: identification 


1. Introduction 


Polycyclic Aromatic Hydrocarbons (PAHs) are 
a family of molecules consisting of carbon and 
hydrogen atoms combined into fused benzenoid 
rings. From a chemical and physical point of view 
they have properties that have led to exciting ap- 
plications in novel materials (Wan et al. 2012; Sul- 
livan & Jones 2008), but at the same time also to 
quite a cautious use on account of their health- 
related impact (Kim et al. 2013; Boffetta et al. 
1997). For astrophysics they play a particularly 
important role since PAHs have been proposed as 
main candidates for carriers of the so-called aro- 
matic infrared bands (AIBs), a series of infrared 
emission features that are ubiquitously observed 
across a wide variety of interstellar objects. These 
emission features are thought to be nonthermal in 
nature and arising from radiative cooling of iso- 
lated PAHs that have been excited by UV radia- 
tion (Sellgren 1984). 


Since they offer such a powerful probe for car- 
bon evolution in space, these bands have been sub- 
ject to extensive experimental and theoretical re- 
search with the ultimate aim being a rigorous iden- 
tification of the molecular structure of the AIB 
carriers. Significant progress has been made in 
this respect with infrared (IR) studies on PAH 
species deposited in a cold (10 K) rare-gas ma- 
trix (for example, Hudgins & Sandford (1998a,b, 


2000)). Although such cooling conditions allow for 
an increase in spectral resolution as compared to 
room-temperature experiments, they lead at the 
same time to matrix-induced effects that are not 
well understood and hard to predict. Gas-phase 
studies are much preferred but have so far predom- 
inantly been restricted to IR absorption studies of 
hot (1000 K) vaporized PAHs (Joblin et al. 1995, 
1994) or at best under room-temperature condi- 
tions (e.g., Pirali et al. (2009); Cané et al. (1997)). 
Due to their low volatility, high-resolution studies 
of low-temperature, isolated PAH molecules has 
for a long time remained out of reach with as a no- 
table exception the cavity ring down spectroscopy 
(CRDS) studies of Huneycutt et al. (2004) on 
small PAHs, although contaminations originating 
from isotopologues and other PAH species or im- 
purities remained a point of concern. 

Recently, we have applied IR-UV double res- 
onance laser spectroscopic techniques on PAHs 
seeded in supersonic molecular beams. In com- 
bination with mass-resolved ion detection these 
techniques allow for recording of mass- and 
conformation-selected IR absorption spectra with 
resonance band widths down to 1 cm~! (Malt- 
seva et al. 2015). Under such high-resolution con- 
ditions, IR absorption spectra of PAHs in the 
3-um region turn out to display an unexpected 
large number of strong bands, and certainly many 


more than expected on the basis of a simple har- 
monic vibrational analysis. Such a conclusion is 
more pertinent as theoretical studies of IR spec- 
tra of PAHs are typically performed at the Den- 
sity Functional Theory (DFT) level, using the 
harmonic approximation for vibrational frequen- 
cies and double harmonic approximation for in- 
tensities, neglecting the effects of anharmonicity. 
In our previous papers we demonstrated that a 
proper treatment of anharmonicity and Fermi res- 
onances indeed leads to predicted spectra that are 
in near-quantitative agreement with the experi- 
mental spectra, both with respect to the frequen- 
cies of vibrational bands and their intensities. 


The shape of the 3-u4m band recorded in as- 
tronomical observations has been found to vary 
within the same astronomical object and between 
different astronomical objects. To account for 
these differences, several explanations have been 
put forward (Sellgren et al. 1990; Tokunaga et al. 
1991; van Diedenhoven et al. 2004). One of the 
suggestions for classifying the shape of this band 
is to interpret it as being associated with emission 
from two components (Song et al. 2003; Candian 
et al. 2012) in which there are contributions from 
different groups of carriers at 3.28 and 3.30 ym. 
Song et al. (2007) propose that these two compo- 
nents originate from groups of PAHs with differ- 
ent sizes, finding support for this in the labora- 
tory high-temperature gas-phase studies of pyrene 
(CieHi9) and ovalene (C32Hi4) for which a blue 
shift of the 3-~m band was observed upon increas- 
ing the size of the PAH (Joblin et al. 1995). An- 
other factor that has been suggested to contribute 
to the apparent two-component appearance of the 
emission are differences in the peripheral struc- 
ture of different PAHs, in particular steric effects 
as occurring for hydrogen atoms at so-called bay- 
sites (Candian et al. 2012). The influence of the 
edge structure has been investigated by means 
of harmonic DFT calculations for large species 
(Bauschlicher et al. 2009, 2008) but systematic ex- 
perimental high-resolution studies are notoriously 
lacking. Similarly, it has been found (Geballe et al. 
1989) that next to the prominent emission band 
at 3.29 um, a broad plateau is present that spans 
the 3.1-3.7 um region and has been indicated as 
the 3-uwm plateau. It has been speculated that 
this plateau might in part derive from anharmonic 
couplings to vibrational combination levels (Alla- 


mandola et al. 1989). However, to what extent 
this explanation can account for the appearance 
of the entire plateau is still far from clear. 


Our previous study aimed at recording spec- 
tra under the highest resolution conditions possi- 
ble and applying the appropriate theoretical treat- 
ment including anharmonic effects and resonances. 
For that reason, we focused on the spectra of 
the linear PAHs naphthalene (CioHg), anthracene 
(Ci4Hi9), and tetracene (CigHj2). As discussed 
above, astronomical spectra likely comprise the 
contributions of a much larger variety of PAHs. 
To advance the interpretation and characteriza- 
tion of these data, we therefore extend our experi- 
mental and theoretical studies to a wider variety of 
condensed and irregular isomers containing up to 
five rings (phenanthrene C14Hj0, benz[ajantracene 
CisHie, chrysene CisHie, triphenylene CisHie, 
pyrene C,gHj9 and perylene C29Hj2). The goals 
of these studies are twofold. Firstly, we aim to 
understand how the effects of anharmonicity ob- 
served for the linear PAHs are affected by geo- 
metrical structure and how this in turn affects the 
appearance of the IR absorption spectra in the 3- 
pum region. Secondly, we aim to uncover general 
trends in band shapes that could provide spectral 
signatures that would allow for a much more de- 
tailed description of the contribution of different 
PAHs. 


2. Methods 


2.1. Experimental techniques 


IR spectra of cold and isolated molecules were 
obtained in molecular beam setup described in 
Smolarek et al. (2011). In these experiments the 
sample of interest was placed in an oven attached 
to a pulsed valve (General Valve). Sufficient va- 
por pressure was found to be obtained when the 
sample was heated to temperatures slightly higher 
than its melting point. Subsequent pulsed expan- 
sion with 2 bars of argon as carrier gas and a typ- 
ical opening time of 200 ys then led to supersonic 
cooling of the sample molecules. 

The ground-state vibrational manifold of PAHs 
was probed by UV-IR ion dip spectroscopy. To 
this purpose two-color resonance enhanced mul- 
tiphoton ionization followed by mass-selected ion 
detection was used to generate an ion signal. In 
these experiments PAHs were electronically ex- 


cited by fixing a frequency-doubled pulsed dye 
laser (Sirah Cobra Stretch) on a) the S; (141) 
+ Sp (1A) 0 - 0 transition at 29329.8 cm~! for 
phenanthrene (Amirav et al. 1984); b) the S, (‘A’) 
+ So (A’) 0 - 0 transition at 26534.4 cm~! for 
benz[ajanthracene (Wick et al. 1993); c) the S; 
(‘Bu) <— So ('Ag) 0 - 0 transition at 28194.3 
cm! for chrysene (Zhang et al. 2012); d) the $; 
(Bou) < So ('Ag) 0 - 0 transition at 27210.9 
cm! for pyrene (which is close to but not ex- 
actly the same as the previously reported value by 
Ohta et al. (1987); e) the S$, (Ai) — So (144) 
transition to the 249 (29867) cm~! (le) level in 
the excited state of triphenylene (Harthcock et al. 
2014; Kokkin et al. 2007); and f) the S; (*B3,) 
< So (1A,) 0 - 0 transition at 24067.4 cm! for 
perylene (de Vries et al. 1991). Molecules in elec- 
tronically excited states were ionized by an ArF 
excimer laser (Neweks PSX-501) in temporal over- 
lap with the electronic excitation laser. 

Prior to the pump and ionization laser beams 
needed to generate this signal, an IR laser beam 
with a linewidth of 0.07 cm7! and a typical pulse 
energy 1 mJ was introduced with a time delay of 
200 ns. The 3 wm beam was generated by differ- 
ence frequency mixing of the fundamental output 
of a dye laser (Sirah Precision Scan with LDS798 
dye) and the 1064 nm fundamental of Nd:YAG 
laser (Spectra Physics Lab 190) in a LiNbO3 crys- 
tal. Resonant excitation of vibrational levels was 
observed by a decrease of the ion signal due to 
depletion of the ground state population, allowing 
for recording of IR absorption spectra. Such IR 
spectra were recorded between 3.17 and 3.40 wm 
(2950 and 3150 cm~'). With the present S/N ra- 
tios, no other IR bands could be observed outside 
this range. 


2.2. Computational methods 


Two types of theoretical IR spectra were pro- 
duced: the standard harmonic vibrational ap- 
proach implemented in Gaussian09 (Frisch et al. 
2009), and an anharmonic vibrational approach 
that employs both Gaussian09 and a locally mod- 
ified version of the program SPECTRO (Mackie 
et al. 2015a; Gaw et al. 1996), referred to in 
this work as G09-h and SP16 calculations re- 
spectively. Both G09-h and SP16 calculations 
start from DFT calculations that employ a sim- 
ilar integration grid as in Boese & Martin (2004), 


the B9-71 functional (Hamprecht et al. 1998), 
and the TZ2P basis set (Dunning 1971), all of 
which has been found to give the best perfor- 
mance on organic molecules (Boese & Martin 
2004; Cané et al. 2007). The SP16 calculations 
utilize the quadratic, cubic, and quartic force 
constants calculated by Gaussian09 and trans- 
formed into Cartesian derivatives as the input for 
SPECTRO, which then performs its own vibra- 
tional second-order perturbation (VPT2) analysis 
(Mackie et al. 2015b,a). SP16 treats the polyads 
of multiple simultaneous resonances (modes falling 
within 200 cm! of each other) and allows for 
the redistribution of intensity among the resonant 
modes, which is a function not currently employed 
in Gaussian09. A more thorough account of the 
theoretical aspect of this work is given in a sep- 
arate publication (Mackie et al. 2016). Harmonic 
calculations in this work were scaled in order to be 
compared with experimental data using a scaling 
factor (sf) of 0.961, while SP16 output does not 
require scaling and can be compared directly. 


3. Results 


In the following we will discuss IR-UV ion dip 
spectra of phenanthrene, benz[ajantracene, chry- 
sene, triphenylene, pyrene and perylene. In the 
first instance, we will present these spectra and 
compare them where possible with previous stud- 
ies, highlighting only salient features and differ- 
ences between the various PAHs. We will make 
extensive comparisons with the results of theo- 
retically predicted spectra. A more detailed dis- 
cussion of these calculations and the assignment 
of each of the observed bands in the experimen- 
tal spectra will be presented in a separate study 
(Mackie et al. 2016). Subsequently, we will make 
a global comparison between the spectra, and dis- 
cuss their implications for extrapolating to the 3- 
pum region of larger PAHs. Finally, we will con- 
sider the current interpretation of astronomical 
PAH data in the light of the present results. 


3.1. Phenanthrene 


The 3-~m absorption spectrum of phenan- 
threne is shown in Figure 1, bottom panel. In total 
twenty-three experimental bands with a linewidth 
2.5-4.5 cm! are observed. The positions and 
strengths of these bands are reported in Table 1. 


The most active region (where the bands are at 
least 40% intensity of the dominant band) is lo- 
cated between 3030 and 3090 cm~! with the most 
intense transition at 3065.3 cm~‘. This region 
is accompanied by two weak bands in the low- 
energy region and five overlapping weak bands in 
the high-energy wing. Table 1 also provides the 
values of bands reported previously in the CRDS 
studies of Huneycutt et al. (2004). Due to better 
cooling conditions in the present experiments our 
absorption spectrum shows bands that are more 
narrow. As a result, nine additional bands can 
be discerned, although we do not find evidence 
for the weak band previously reported at 3025.3 
cm~+. Overall, the agreement between the two 
sets of data is good with the largest deviation for 


most of the bands not exceeding 2 cm™. 


Phenanthrene belongs to the Cz, point group. 
Within the harmonic approximation (Fig. 1, top 
panel) only six modes are therefore IR active in 
the 3-um region. The experimental spectrum, on 
the other hand, displays twenty-three bands and 
is a clear indicator of the failure of the harmonic 
approximation. We conclude that, similar to lin- 
ear PAHs (Maltseva et al. 2015), the CH-stretch 
region is dominated by Fermi resonances. Pre- 
viously, we showed (Mackie e¢ al. 2015b) that in 
order to obtain good agreement between exper- 
imental and predicted spectra, VPT2 treatments 
are required that take intensity sharing and polyad 
resonances as implemented in SP16 into account. 
With such calculations, the redistribution of inten- 
sities over fundamental and combination bands — 
in this case overtones are not accessible due to 
symmetry restrictions — of the same symmetry 
leads to observable activity of many additional 
bands. Comparison between the 3-~m absorp- 
tion spectrum of phenanthrene predicted by SP16 
and convolved with a 1 cm~! Gaussian line shape 
(Fig.1, middle panel) with the experimentally ob- 
tained spectrum indeed shows good agreement. 


3.2. Pyrene 


Figure 2 (bottom panel) displays the IR absorp- 
tion spectrum of pyrene in the 2950-3150 cm~+ 
region as measured in the present study together 
with spectra predicted by G09-h and SP16 cal- 
culations. Compared to phenanthrene, the dom- 
inant vibrational activity is observed in a more 
compact region of 3045-3065 cm~! with the maxi- 


mum intensity at 3049.8 cm~!. The experimental 
spectrum shows fourteen relatively narrow bands 
with a linewidth of 1.7-3.5 cm7! whose positions 
are given in table 1. This table also contains the 
positions of the ten bands reported previously by 
Huneycutt et al. (2004), which agree well with 
the present observations (deviations less than 1 
cm~'). Our S/N ratio, on the other hand, allows 
us to identify weak bands at 3044, 3087.8, and 
3096 cm! that could previously not be discerned, 
but does not confirm the previously reported band 
at 3083.9 cm~!. Due to the Do, symmetry of 
pyrene, only five IR active CH-stretch modes are 
expected in the harmonic approximation (Fig.2, 
top panel), which is clearly at odds with the ex- 
periment. Indeed, the polyad VPT2 anharmonic 
calculations (Fig.2, middle) are required to bring 
experiment and theory into better agreement. 


3.3. Benz/ajanthracene 


The IR absorption spectrum benz[ajanthracene 
in the 2950-3150 cm~* region is shown in Figure 3 
together with the spectra predicted by harmonic 
and anharmonic calculations. The main spectrum 
feature covers quite a large frequency range and 
lacks the narrow bands observed for the other 
PAHs reported here. Instead a broad structure 
between 3030 to 3090 cm~! with the highest in- 
tensity at 3063.8 cm7! is seen. In this structure 
only five bands can be distinguished clearly due 
to the overlap of a large number of blended tran- 
sitions in the present resolution. 

The presently obtained spectrum is the first 
high-resolution gas-phase IR spectrum that has 
been reported for benz[ajanthracene. We there- 
fore compare in Table 1 the line positions as de- 
rived from Figure 3 with the Matrix Isolation 
Spectroscopy (MIS) data of Hudgins & Sandford 
(1998b). We find larger deviations of up to 3cm7! 
which likely find their origin in matrix-induced 
perturbations, difficulties in determination of the 
center of bands due to their shape (as occurs for 
the most intense band), and the noise level (as oc- 
curs for the weak bands). Interestingly, our spec- 
trum does not give an indication for the presence 
of a band at 3119.5 cm~! that was reported in the 
MIS studies, possibly originating from the lack of 
mass selection in the MIS experiments. 

The large frequency range over which vibra- 
tional activity is observed is in line with the rel- 


atively low symmetry of the molecule (C,). As a 
result, all twelve CH-stretch modes are formally IR 
active, although the harmonic calculations predict 
that only six of them have appreciable intensity 
(Fig.3, top). The SP16 calculations, on the other 
hand, match the experimental spectrum very well. 
They predict activity of a plethora of vibrational 
transitions as is clear when the stick spectrum un- 
derlying the green trace in Fig. 3 (middle panel) 
is inspected, and thereby confirm the conclusions 
that the the broad appearance of the spectrum is 
due to the overlapping of many anharmonic bands. 


3.4. Chrysene 


Chrysene (CigHj2) is a highly symmetric iso- 
mer of benz[ajanthracene (C2), versus C,). Its ex- 
perimentally obtained and theoretically predicted 
IR absorption spectra are shown in Figure 4. In 
contrast to benzlajanthracene, the experimental 
spectrum displays a larger number of resolved 
bands that have line widths which range from 3-7 
cm~'. Remarkably, major vibrational activity is 
found in the 3015-3115 cm~! range, which is even 
larger than for benz[ajanthracene, with the most 
intense band at 3063 cm7!. In this respect, it is 
interesting to note the presence of bands above 
3090 cm~! where no activity was observed for 
benz/ajanthracene. 

Table 2 reports the line positions of the fif- 
teen bands that can be distinguished and com- 
pares them with previously reported MIS positions 
(Hudgins & Sandford 1998b). This comparison 
shows deviations less than 2 cm~! caused mainly 
by the inability to accurately identify line posi- 
tions of close lying bands under the relatively low- 
resolution conditions of the MIS experiment. The 
gas-phase data furthermore identify eight previ- 
ously unreported transitions, but does not show 
evidence for the MIS band at 3134.9 cm7!. Sym- 
metry considerations lead to the conclusion that 
at most six CH-stretch modes could be IR-active 
(Fig.4, top panel) and that they are localized in a 
relatively narrow frequency range. SP16 calcula- 
tions (Fig.4, middle panel) once again emphasize 
the important role of anharmonicity and Fermi 
resonances, and give rise to a predicted spectrum 
that is in good agreement with the experiment. 


3.5. Triphenylene 


Triphenylene completes the series of PAHs with 
four aromatic rings. For this molecule no high- 
resolution gas-phase vibrational spectra have been 
reported previously. The 3-zm band recorded here 
using UV-IR ion dip spectroscopy is depicted in 
Figure 5 (bottom panel). Activity is observed over 
a relatively small range of 3030-3105 cm~! with 
three strong features at 3040, 3075 and 3100 cm7! 
and the strongest band at 3101.8 cm~!. Above 
3105 cm7!, one band at 3143.2 cm~! appears to 
be present, there is no indication for bands be- 
low 3025 cm7!. Overall, sixteen bands with an 
irregular shape and a minimum linewidth of 2.9 
cm! can be resolved. Considering their profile, 
it is likely that most of these bands actually con- 
sist of several overlapping transitions. Previous 
MIS studies (Hudgins & Sandford 1998b) compare 
well with the present data (see Table 2) with dif- 
ferences in band positions not exceeding 2 cm™!. 
In the present study significantly more bands are 
observed than in the MIS studies, and the MIS 
3116.9 cm~! band is not confirmed. 


Triphenylene belongs to the D3, symmetry 
point group and because of this high symmetry 
only four pairs of doubly degenerate CH-stretch 
modes are IR active in the harmonic approxi- 
mation. As Gaussian 09 can not perform the 
anharmonic analysis of molecules with D3, sym- 
metry, the calculation was done on triphenylene 
with a reduced C2, symmetry caused by a small 
perturbation of the masses of two opposing car- 
bon atoms, from 12 to 12.01 au (see Mackie et al. 
(2016) for further details). 

Under such conditions a spectrum is predicted 
(Fig. 5, top) with only two strong bands, which is 
clearly at odds with the experiment. SP16 calcu- 
lations incorporating anharmonicity, on the other 
hand, redistribute the intensities among funda- 
mentals and combination bands with equal sym- 
metries and lead to a predicted absorption spec- 
trum that resembles the experimentally observed 
one (Fig.5, middle panel). Nevertheless, a detailed 
comparison of the two spectra does lead to the 
conclusion that in the regions above 3105 and be- 
low 3025 cm~! quite a larger activity is predicted 
by the calculations than observed experimentally. 
In fact, similar observations can also be made 
for chrysene (regions above 3110 and below 3020 


cm~') and to some extent for benz[a]anthracene 
(region above 3090 cm~') although in these cases 
the differences are not as obvious as in the case of 
triphenylene. 


3.6. Perylene 


Perylene is the most condensed five-aromatic 
ring PAH system. The experimental IR absorp- 
tion spectrum of this molecule in the 3-yum re- 
gion is shown in the bottom panel of Figure 6 to- 
gether with predicted spectra at the GO09-h level 
(top panel). The spectrum is quite compact and 
dominated by a group of bands located in 3056- 
3070 cm—! with the most intense transition oc- 
curring at 3063.8 cm~!. At the high-energy side 
a weak, sharp band is observed at 3118 cm! and 
a broad structure that appears to consist of three 
overlapping bands at 3092.8, 3095.9, 3098.2 cm7!. 
Overall, nine bands with a linewidth of >3 cm7! 
are reported in Table 2. Perylene has previously 
been studied with CRDS (Huneycutt e¢ al. 2004), 
but the spectrum reported at that time clearly has 
a lower S/N ratio. This is most likely the reason 
that a comparison between the two spectra is not 
as favorable as for the other molecules. For exam- 
ple, we do not find indications for the previously 
reported bands at 3044.9 and 3070.4 cm7! and 
cannot confirm the band at 3022 cm7!. More- 
over, our cooling conditions enable us to conclude 
that the bands reported at 3096.1 cm~! and 3065.3 
cm~* (CRDS) consist of three and two bands, re- 
spectively. Perylene possesses D2, symmetry and 
according to the harmonic approximation only 6 
CH-stretch modes are IR active. The other ob- 
served bands originate from anharmonic activities 
—i1.e., resonances. 

Unfortunately, it was not possible to perform 
an anharmonic analysis on perylene as Gaussian09 
was unable to produce reasonable cubic and quar- 
tic force constants. Those unreasonable force con- 
stants lead to a significant overestimation of the 
anharmonic corrections of band positions, with 
many bands predicted to occur at higher energies, 
which is contrary to what is expected. It would 
appear that these problems depend on the level of 
theory that is used since calculations on perylene 
using the B3LYP functional with a 4-31G basis 
set lead to cubic and quartic force constants that 
seem much more reasonable. However, this level 
of theory is too low to allow for meaningful com- 


parisons with the experimental spectra. Further 
studies on the functional/basis set dependence of 
the quartic force field of large molecules such as 
perylene would thus be of interest, but fall out- 
side the scope of the present studies. 


4. Discussion 


The high-resolution IR absorption spectra of 
the presently studied set of PAHs together with 
those of the linear acenes studied previously 
(Maltseva et al. 2015) provide a comprehensive 
view on key factors that determine the appear- 
ance of the 3-4m band in PAHs. The number 
of CH-oscillators obviously controls the number 
of normal modes, and this clearly implies that 
compact species always have less active modes in 
the 3-um region than more extended molecules of 
the same symmetry and with the same number 
of rings. However, from the comparison of 3-~m 
bands of PAHs with different molecular structure 
it becomes clear that the frequency region over 
which activity is observed and the distribution of 
intensity over this region is not as dependent on 
the number of CH-oscillators as one might expect. 


For example, four-ring pyrene consists of six- 
teen carbons and ten hydrogens and its 3-y~m 
bands looks quite similar to anthracene (Fig.7), 
another molecule with ten CH-oscillators and Do, 
symmetry. Both 3-u~m bands are compact, the 
most active region occurring in a range of 22 
and 30 cm™!, respectively. However, a notice- 
able blue-shift of the strongest transition of an- 
thracene (3071.9 cm~!) with respect to pyrene 
(3049.8 cm~') is observed. Phenathrene is an iso- 
mer of anthracene, but despite just a small change 
in molecular structure, its 3-um band is very 
different (Fig.7). The most active region spans 
60 cm~! with the most intense band at 3065.3 
cm‘. These differences likely originate from the 
nonlinear structure of phenanthrene. Significant 
changes occur in the periphery of the molecule. 
Anthracene has two solo-hydrogens and two quar- 
tets. Phenanthrene, on the other hand, has two 
duo hydrogens instead of two solos, and also two 
quartets, but conversely the quartets each have a 
hydrogen in the bay-region (Fig. 8). As a result, 
steric effects increase the vibrational frequencies 
of these CH-oscillators (Bauschlicher et al. 2009). 
Interestingly, Figure 9 shows that the 3-wm bands 


of phenanthrene resembles to a large extent that 
of benz[ajanthracene. The most intense transi- 
tions of phenanthrene and benz|ajanthracene are 
at 3065.3 and 3063.8 cm~!, the most active regions 
are 3030-3090 and 3035-3097 cm~', respectively. 
This can be rationalized on the basis of their sim- 
ilar peripheral structure, the only difference be- 
tween the two molecules being the additional two 
solo hydrogens in benz[ajanthracene, which do not 
contribute significantly to the spectrum. 


The 3-um absorption of pyrene looks very sim- 
ilar to the absorption of the five-ring PAH with 
twelve hydrogens - perylene (Fig.9). The most ac- 
tive modes of perylene are concentrated in a nar- 
row frequency region of 15 cm~!. This relatively 
small range follows quite nicely from the molecu- 
lar and, in particular, the peripheral structure of 
perylene. All hydrogens in perylene come as trios 
and thus contribute to IR absorption at similar 
frequencies. Perylene and pyrene indeed have a 
similar band at 3064 cm~! which is attributed to 
modes localized in both cases in trio hydrogens. 
Pyrene does not possess bay-hydrogens, but pery- 
lene does. As confirmed by the calculations this 
explains the differences in activity in the high- 
frequency regions of the spectra of the two com- 
pounds. 


From the subsequent discussion (vide infra) it 
appears that molecules with the same type of hy- 
drogens (solos, duos, trios, or quartets) display 
activity in a similar frequency range. Moreover, 
from the comparison of our experiment with DFT 
theory, we suggest that the frequencies of normal 
modes that primarily involve motion of the same 
type of hydrogens are restricted to rather narrow 
frequency regions with frequencies that increase in 
the order of solos, duos, trios, to quartet (Mackie 
et al. 2016). 


Another effect induced by steric hindrance asso- 
ciated with bay regions is best illustrated by com- 
paring the 3-~m absorption of tetracene with its 
isomers in Fig.10. Benz[ajanthracene has a pair 
of hydrogens (one is from a quartet and another 
one is a solo hydrogen) in a bay-site. Its 3-~m 
band has a maximum at 3063.8 cm~!, close to the 
most intense transition of tetracene (3061.1 cm~') 
and major activity over 3035-3097 cm~!. While 
the strongest transition of chrysene (3063 cm~!) is 
not affected by the structural changes and remains 
very close to tetracene, the high-energy part of 


the spectrum undergoes changes. Chrysene shows 
prominent high-energy features in the 3090-3120 
cm~+ region where tetracene has considerably less 
activity (Fig.10). We conclude that this activity 
should be attributed to bands involving modes 
with considerable contributions of hydrogens in 
the bay region as these have intrinsically higher vi- 
brational frequencies because of steric hindrance. 
The conclusion is supported by the SP16 calcula- 
tions. It is found computationally that all reso- 
nances in this region involve the CH-stretch mode 
of chrysene located in a bay region. All hydrogens 
of triphenylene are incorporated into three quar- 
tets with half of them in bay-sites. In contrast 
to the other PAHs with quartets, triphenylene has 
the most intense band at 3101.8 cm~! while for 
the other compounds it is found between 3063- 
3066 cm~! (Fig.10). This is a clear demonstration 
of the presence of bay-site hydrogens. The band 
at 3076.2 cm~! most likely dominantly involves 
asymmetrically stretching of bay-hydrogens. The 
3030-3060 cm~! region can nicely be explained by 
anharmonic activity (Mackie et al. 2016). Such a 
comparison clearly shows that vibrations involving 
hydrogens in bay-site are blue-shifted. 


The symmetry of the molecule controls the 
number of [R-allowed transitions. These transi- 
tions are the source of intensity for combination 
bands of the same symmetry that normally would 
be ’dark’ but can acquire intensity through Fermi 
resonance. When comparing isomers with differ- 
ent symmetries, for example anthracene and phe- 
hanthrene, the isomer with a higher symmetry 
(anthracene) is thus expected to show a smaller 
number of fundamental and combination bands as 
is indeed confirmed by our experiments (Fig.7). 
Therefore, molecules with low degree of symmetry, 
like benz[aJanthracene, demonstrate broad absorp- 
tion because of the overlap of anharmonic bands 
(Fig.3). 

Another important conclusion that can be 
drawn from the present experiments concerns the 
frequency range over which anharmonicity is ac- 
tive in the 3-4m region. We have acquired high- 
resolution IR absorption spectra for PAHs with 
two to five rings in all possible geometrical ar- 
rangements, but for all these species vibrational 
activity quickly disappears below, roughly speak- 
ing, 3000 cm~! (above 3.33 um). This conclusion 
is supported by our SP16 calculations that also 


show very limited activity in this region. More- 
over, in the SP16 calculations an important pa- 
rameter that controls the number of states to be 
included in resonanaces/polyads is the maximum 
energy separation of the states. In all our cal- 
culations we find that a maximum separation of 
200 cm~! is more than adequate. On the basis 
of the scaled harmonic calculations in which the 
lowest frequencies occur around 3025 cm! one 
can then conclude that the maximum wavelength 
range over which combination band activity might 
be observed extends only up to 2825 cm~! (3.54 
pm), and this number is most likely even an over- 
estimate. Increasing the size of a PAH will thus 
increase the density of ‘dark’ states that might po- 
tentially couple to the ‘bright’ zeroth-order states, 
but the present study strongly suggests that this 
increased density has a very limited effect on the 
range of the region over which vibrational activity 
is observed. 


5. Astrophysical implications 


The present study confirms the dominant role 
of anharmonicity in determining the shape and 
strength of the features in the CH-stretching re- 
gion. Predictions by harmonic calculations thus 
are far from adequate for obtaining a fundamen- 
tal understanding of how the 3-ym region reflects 
the chemical composition and evolution of astro- 
nomical objects, and for trying to understand sec- 
ondary features around the 3-um feature. They 
also indicate that in principle there are many more 
leads available for identifying single PAH species 
as intensity is distributed over many more transi- 
tions than only the fundamental CH-stretch tran- 
sitions. 

One of the secondary features that has been dis- 
cussed extensively in literature concerns the two- 
component emission character of the 3-u4m band 
(Song et al. 2003). Our results clearly show that 
bay-hydrogens induce intensity at the high-energy 
side, and thus fully support the previous sugges- 
tion that the 3.28-~4m emission band is associ- 
ated with modes involving bay-hydrogens as oc- 
curring in more extended PAHs, while the 3.30-um 
emission band derives from more compact PAH 
structures (Candian et al. 2012). The size of the 
PAH, which also has been put forward as the pri- 
mary cause for the two-component emission (Song 


et al. 2007), appears to be only of secondary influ- 
ence. Nevertheless, our results also unmistakably 
demonstrate that activity in the high-frequency 
range is not only induced by bay-hydrogens, but 
also occurs for PAHs without bay-sites. The finer 
details of this part of the spectrum thus are the re- 
sult from a subtle interplay between the effects of 
steric hindrance as encountered for the bay hydro- 
gens and anharmonicities. This implies that one 
should be cautious in directly relating the ratio of 
the intensities of the two bands to abundances of 
compact and extended PAH structures. 


Another secondary feature that has attracted 
considerable attention is the 3-um plateau span- 
ning the 3.1-3.7 wm (2700-3200 cm~') region (Al- 
lamandola et al. 1989; Geballe et al. 1989). Our 
spectra show that the intensity in the high-energy 
part of this plateau is derived from Fermi res- 
onances between fundamental CH-stretch transi- 
tions and the plethora of combination bands that 
are present in this region. However, such an ex- 
planation, which previously has also been put for- 
ward as a possible cause for the lower-energy side 
of the plateau (Allamandola et al. 1989), is not 
supported by the present results which show only 
very limited activity in the region above 3.3 um 
(below 3000 cm~'). Above, we have reasoned that 
also for larger PAHs we do not expect an increased 
activity in this region. We therefore conclude that 
in this region IR absorption has another origin. 

A final feature of interest is the 3.40-~m fea- 
ture on the 3-um plateau. Several explanations 
have been put forward to account for its presence 
ranging from hot bands of aromatic CH-stretch 
transition v = 2 > 1 shifted due to anharmonic 
effects (Barker et al. 1987) to CH-stretch modes 
in methylated (Joblin et al. 1996; Pauzat et al. 
1999) and superhydrogenated (Bernstein et al. 
1996; Steglich et al. 2013; Sandford et al. 2013; 
Wagner et al. 2000) PAHs , but none of these have 
so far been confirmed. The overall picture that 
emanates from our high-resolution studies, involv- 
ing in particular the role of anharmonicity, sug- 
gests that an explanation based on hydrogenated 
and alkylated PAHs is attractive. Methylated and 
hydrogenated PAHs have been known to be IR 
active at much lower frequencies than the fully 
aromatic systems and might thus contribute to 
the 3-um plateau. Assuming that anharmonicity 
plays a similar role in such alkylated PAHs, this 


would then also explain the extended range of the 
3-m plateau which cannot be explained merely 
in terms of bare PAHs. To find further support 
for such a conclusion we are presently performing 
high-resolution IR absorption studies on an exten- 
sive series of hydrogenated and alkylated PAHs. 


6. Conclusions 


In this work we have presented molecular beam 
IR absorption spectra of six condensed PAHs in 
the 3-ym region using IR-UV ion dip spectroscopy. 
The present results and the results on linear acenes 
(Maltseva et al. 2015; Mackie et al. 2015b) show 
that anharmonicity indeed rules the 3-wm region, 
with the fraction of intensity not associated with 
fundamental transitions easily exceeding 50%. 
Anharmonicity-induced transitions are more the 
rule than the exception and should explicitly be 
taken into account in the interpretation of astro- 
nomical data in the 3-4m region. A proper incor- 
poration of resonances has been shown to yield 
predicted spectra that are in semi-quantitative 
agreement with the experiments. Such calcula- 
tions may therefore lead the way for furthering 
our understanding on the influence of larger PAHs 
which are not amenable to similar experimental 
high-resolution studies. 


Our work shows that the observed abundance 
of combination bands is mostly concentrated in 
the low-energy part (< 3100 cm~') of the CH- 
stretch region and originates from combinations of 
CC-stretch and CH in-plane bending modes. This 
anharmonic activity can be partly responsible for 
the 3-ym plateau observed by astronomers al- 
though both experiment and theory put into ques- 
tion a scenario in which the plateau is solely at- 
tributed to Fermi resonances between fundamental 
modes and such combination bands. In this re- 
spect, a more important role than assumed so far 
of hydrogenated and alkylated PAHs in combina- 
tion with anharmonic effects appears to provide a 
highly-interesting alternative to pursue. Activity 
at the high-energy part of the 3-wm band has been 
demonstrated to derive from the presence of bay- 
hydrogens, but not solely as also anharmonicity 
in PAHs without bay-hydrogens induces activity 
in this region. Further studies on larger compact 
PAHs should provide further means to distinguish 
between the relative importance of the two effects. 
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Such studies are presently underway. 


Our studies show that the vibrational activity 
in the CH-stretch fingerprint region is strongly 
linked to details of the molecular structure. It 
is therefore not only different for molecules with 
different chemical structures, but also for different 
isomers. We have demonstrated that the periph- 
eral structure of the molecules plays a dominant 
role in the appearance of the 3-m band. Studies 
in which the 3-ym region is correlated with the 
*periphery-sensitive” 9-15 jum region are thus of 
significant interest and could further elucidate the 
exact composition of the carriers of these bands. 
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Table 1: Frequencies (cm~') and intensities of phenanthrene, pyrene and benz[ajantracene obtained in our 
study and compared to the data from cavity ring down spectroscopy (CRDS)(Huneycutt et al. 2004) and 
matrix-isolation spectroscopy (MIS) data (Hudgins & Sandford 1998b). 


Phenanthrene (C14Hi0) Pyrene (Ci6Hio) benz{alantracenee (CigHj2) 
this work CRDS this work CRDS this work MIS 

freq. rel. int. freq. freq. rel. int. freq. freq. rel. int. freq. 

3015.5 0.25 3044 0.34 3045.9 3017.1 0.2 3016.5 

3018.6 0.26 3018.7 3049.8 1 3051.7 3036.8 0.57 

3032.8 0.68 3034.7 3052.9 0.6 3047.8 0.73 3044.5 

3037.6 0.47 3038.5 3055.6 0.29 3063.8 1 3064.1 

3042.3 0.45 3057.5 0.13 3079 0.69 3077.6 

3042.4 0.45 3044.3 3059.5 0.39 3058.1 3087.8 0.37 

3047.9 0.55 3063.3 0.56 3061.3 

3050.3 0.4 3049.1 3064.9 0.59 3064.5 

3056.4 0.88 3057.3 3067.1 0.12 3067.3 

3057 0.9 3058.5 3071.9 0.11 3071.7 

3061.8 0.5 3087.8 0.08 3090 

3065.3 1 3066.9 3096 0.12 3094.8 

3069 0.57 3108.9 0.1 3108.4 

3071.6 0.64 3072.3 3118.7 0.12 

3075.4 0.89 3076.7 

3081.5 0.5 

3082.8 0.53 

3084.3 0.55 3083.9 

3091.8 0.3 

3094 0.33 3093.7 

3100 0.28 3101.3 

3109.6 0.25 3109.7 

3116.6 0.24 
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Table 2: Frequencies (cm!) and intensities of chrysene, triphenyene and perylene obtained in our study 
and compared to the data from cavity ring down spectroscopy (CRDS)(Huneycutt et al. 2004) and matrix- 
isolation spectroscopy (MIS) data (Hudgins & Sandford 1998b). 


Chrysene (CigHj2) Triphenylene (CigHj2) Perylene (CigHj2) 

this work MIS this work MIS this work CRDS 
freq. rel. int. freq. freq. rel. int. freq. freq. rel. int. freq. 
3010.1 0.2 3035.3 0.64 3048.9 0.3 3050.3 
3020.3 0.5 3021.9 3039.4 0.5 3038.8 3057.4 0.5 
3035.1 0.5 3033.4 3040.6 0.35 3059.5 0.48 3060.9 
3043.5 0.32 3046.4 0.4 3063.8 1 3065.3 
3049 0.45 3050.8 0.57 3049.4 3066.3 0.85 
3053.7 0.75 3054.1 3058.2 0.34 3092.8 0.24 
3058.8 0.56 3065 0.34 3095.9 0.21 3096.1 
3063 1 3063.9 3070.3 0.31 3098.2 0.24 
3070.5 0.71 3073.8 0.62 3074.7 3118 0.3 
3073.8 0.93 3076.2 0.85 
3078.4 0.78 3079.4 3084.3 0.36 
3091 0.23 3088.1 0.35 
3095.7 0.65 3093.8 3091.5 0.43 3090.9 
3102.3 0.4 3098.7 0.68 3098 
3107.4 0.54 3106.8 3101.8 1 

3143.2 0.2 
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Fig. 1.— The IR absorption spectrum of phenanthrene as predicted by G09-h (scaling factor sf= 0.961) and 
SP16 calculations (not scaled) together with the molecular beam gas-phase spectrum as measured in the 
present experiments. Number of hydrogens is mentioned in the following order: solo/duo/trio/quartet /bay 
H’s. See Fig.8 for more details. 
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Fig. 2— The IR absorption spectrum of pyrene as predicted by G09-h and SP16 calculations together with 
the molecular beam gas-phase spectrum as measured in the present experiments. 
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Fig. 3.— The IR absorption spectrum of benz[ajanthracene as predicted by G09-h and SP16 calculations 
together with the molecular beam gas-phase spectrum as measured in the present experiments. 
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Fig. 4.— The IR absorption spectrum of chrysene as predicted by GO9-h and SP16 calculations together 
with the molecular beam gas-phase spectrum as measured in the present experiments. 
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Fig. 5.— The IR absorption spectrum of triphenylene as predicted by G09-h and SP16 calculations together 
with the molecular beam gas-phase spectrum as measured in the present experiments. 
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Fig. 6.—\ The IR absorption spectrum of perylene as predicted by GO09-h calculations together with the 
molecular beam gas-phase spectrum as measured in the present experiments. Calculations at the SP16 level 
lead to incorrect results (see text) and therefore are not reported. 
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Fig. 7.— Experimental IR absorption spectra of pyrene, anthracene and phenanthrene from the bottom to 
the top, respectively. Number of hydrogens is mentioned in the following order: solo/duo/trio/quartet /bay 
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Fig. 8.— Schematic representation of bay-region, solo, duo and quartet hydrogens using examples of an- 
thracene (bottom) phenanthrene (top). 
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Fig. 9.— Experimental IR absorption spectra of pyrene, perylene, benz[a]anthracene and phenanthrene from 
the bottom to the top. 
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Fig. 10.— Experimental IR absorption spectra of tetracene, benz[ajantracene, chrysene, triphenylene from 
the bottom to the top. 
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